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Abstract 

Interaction of a two-level atom with a single mode of electromagnetic 
field including Kerr nonlinearity for the field and intensity-dependent 
atom-field coupling is discussed. The Hamiltonian for the atom-field sys- 
tem is written in terms of the elements of a closed algebra, which has 
SU(l,f ) and Heisenberg-Weyl algebras as limiting cases. Eigenstates and 
eigenvalues of the Hamiltonian are constructed. With the field being in 
a coherent state initially, the dynamical behaviour of atomic-inversion, 
field-statistics and uncertainties in the field quadratures are studied. The 
appearance of nonclassical features during the evolution of the field is 
shown. Further, we explore the overlap of initial and time-evolved field 
states. 
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1 Introduction 



Interaction of a single mode of electromagnetic field with a two-level atom is 
the simplest problem in matter-radiation coupling. A model for the interaction, 
introduced by Jaynes and Cummings [1], treats the atom as a dipole placed 
in an external field. The Jaynes-Cummings (JC) model has provided a lot of 
impetus for theoretical explorations and experimental verifications [2, 3, 4, 5, 6]. 
The Hamiltonian for the model is 

Hjc =Ho+ g{a^a- + da+). (1) 

Here Hq, the Hamiltonian for the atom and field in the absence of interaction, is 
oja^a + \ vaz. The atomic transition frequency is v, the field frequency is a; and 
g is the coupling constant. We denote the creation operator by a\ annihilation 
operator by a, and their action on the basis states of the harmonic oscillator are 



a\n) = ^/n|n-l), (2) 
a|0) = 0, (3) 
a^\n) = Vn+ l\n+ 1). (4) 

The atom has two-levels, \g) and |e), the ground and excited states. The oper- 
ators az , <J+ and cf- act on the states as given below: 

<^A9) = -\9) (5) 
a,|e) = \e) (6) 

<r±\g) = ^|e) (7) 

a±|e) = i^ls) (8) 

(9) 

Atomic inversion, defined as the difference in probabilities for the atom 
to be in the excited and ground states, as predicted by JC model is a sum 
of quasiperiodic functions with incommensurate frequencies. The model pre- 
dicts collapses, revivals and ringing revivals in the time-development of atomic 
invcrsion[7, 8, 9, 10, 11] The revival phenomenon is entirely quantal, and hence 
the model is very important for testing predictions of quantum theory. The 
model has been generalized in very many ways. We list the respective Hamilto- 
nian for some of the models: 
1) Buck-Sukumar model[12] 

Hjc = Hq + g{V a^daa+ + a^Va^aa^) (10) 

With this particular form of intensity-dependent coupling the atomic-inversion 
is a sum of periodic functions. This model is very interesting as it can be writ- 
ten as a combination of generators of SU(1,1) algebra[13]. Generalization to 



2 



two-photon case, where two-photons are absorbed or emitted in atom-field in- 
teractions, has been done. Once again the SU(1,1) algebraic structure in the 
model is used to solve the problem [14]. 

2) Kerr-nonlinearity[15, 16, 17] 

HKerr = Hjc + xa^^a^ (11) 

This is an efi'ective Hamiltonian for a system in which the electromagnetic field 
mode is excited in a Kerr medium. The medium is modelled as an anharmonic 
oscillator [18, 19]. 

3) Dicke-Tavis-Cummings model[20, 21] 

In this model, the interaction between field and a group of two-level atoms 
is considered and the Hamiltonian is 

Hdtc = ujo^a + ^ (Ji^z + Interaction part (12) 

i 

4) Nonlinear Jaynes-Cummings model 

By including the motion of atom in the external field, the coupling is made 
position dependent. This offers enormous possibilities to tailor the form of 
atom-field interaction [22, 23]. The general form for the Hamiltonian is 

Hnl = Ho+ g{f{a^a)br + adjoint) (13) 

Here f{a)a) is an operator- valued function of the number operator a) a and m 

is an integer. 

The objective of the present paper is to study the dynamics of a two-level 
atom interacting with a single mode of electromagnetic field. The interaction is 
governed by the Hamiltonian 

+ xa^'^c? + gu){\/l + ka^aaa^ + a^-\/l + ka^aa-). (14) 

This Hamiltonian is an example for nonlinear JC model including Kerr term. 
Note that we have scaled the coupling constant g hy u). The parameter r is - 
and the coTipling constant for the Kerr term is % = kuj. This will simplify the 
expressions we derive in sequel. In particular, we set oj = 1 which amounts to 
studying a new Hamiltonian The speciality of the Hamiltonian H is that 
it becomes Hjc when the parameter k is set equal to zero. Further, the usual 
Holstein-Primikoff realization is obtained when k = 1 so that the interaction is 
approximately given by g{V a^aaa+ + adjoint), which is same as the interaction 
studied in Ref.[13]. Yet another form of interaction occurs when k <^1. If the 
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photon number distribution is such that kn <^ 1 for aU n under the peak of the 
distribution, then Eq. 14 leads to 

H ^LO [a^a + ka^^a^ + g[{l + a)a/2)a(j+ + adjoint]] . (15) 

Further, when A: <C 1 so that we can neglect kn in comparison to unity but 
retain kv? , we arrive at Hjc with an additional Kerr term. This system has 
been studied in Rcf . [24] . Many well studied systems arc thus special cases of the 
Hamiltonian considered for discussion in the present paper. The organization 
of the paper is as follows. In section II we study the algebraic aspects of the 
Hamiltonian and construct the eigenvectors and eigenvalues. Section III is de- 
voted to study the atomic inversion and approximate expressions are obtained 
for first collapse and revival periods when the field is in a coherent state. The 
time-development of field statistics is discussed in section IV and the results are 
summarized in Section V. 

2 Model Hamiltonian and its properties 

In this section we study the algebraic aspects of the generalized Hamiltonian 
given in Eq. (14). We introduce a set of operators which are closed under com- 
mutation. This algebraic structure is exploited to determine the time-evolution 
operator to evolve the initial state of the atom-field system. The eigenvalues 
and the corresponding eigenvectors are constructed. 

2.1 Eigenvalues and eigenvectors 

The Hamiltonian given in Eq.l4 is written as 

H = ujK+K_ + '^<j,+g{K+a-+K_a+), (16) 

wherein we have set K- = Vl + ka^aa and = a^y/l + ka)a. Further, we 
assume that k is nonnegative and restricted to take values less than or equal 
to unity. Formally, the Hamiltonian H has the same structure as Hjc with 
a and a) replaced by K- and However, the former corresponds to Kerr- 
type medium with intensity-dependent coupling for atom-field interaction. The 
difference is very clear in the commutation relations among the operators. From 
the realisation of K- and Kj^ in terms of and a, we arrive at 

[K^,K+] = 2Ko, 

[Ko.K±] = ±kK±. (17) 

The operator Kq is ka^a+^. Thus, the operators K-,K+ and Kq form a closed 
algebra. It worth noting that the commutation relations define the SU(1,1) 
algebra when k = 1. On the other hand, to get the well-known Heisenberg-Weyl 
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algebra generated by , a and the identity / we set A: = 0. Two different algebras 
are realized depending on the value of k and hence the algebra of K- , K+ and Kq 
is said to be an "interpolating algebra" . An invariant operator, which commutes 
with K± and Kq, for the algebra is given by - {k/2){K_K+ + K+K_). The 
coherent states corresponding to this algebra and their Hilbert space properties 
are known [25]. 

The atom-field evolution is studied in the space of |e,n) and \g,n), where 
n = 0, 1, 2, The state |e, n) means that the atom is in the excited state |e) and 
the field in the excited state \n). The states |e, n) and \g, n) are eigenstates of 
u!K+K- + ^az and the respective eigenvalues are fe,n = {n+kn'^ — kn)u) + ^ and 
£g,n = {n + kv? — kn)u> — f . The Hamiltonian H admits a constant of motion N 
such that the commutator [TV, H] vanishes. Explicitly, N = K+K^ +2KQa+<T^. 
Note that when k = 0, N becomes a^a + a+a-, the constant of motion for Hjc- 

The interaction part of H is such that the state |e, n) is taken to l^, n + 1) 
and vice versa, during the evolution of the atom- field system. Thus, the entire 
Hilbert space is split into subspaces spanned by \e,n) and \g,n + 1) and the 
dynamics confined to individual subspaces. In one such subspace, specified by 
the value of n, the Hamiltonian matrix is 



„ ^ / + I + kn^ - kn)io + A/2 g ^/{T+kn){l + n) \ 

V 5V(l + A;n)(l + n) {n+\+kn'^ + kn)uj - A/2 y " ^ ' 

The detuning parameter A is (r — 1)(jO. The eigenvalues of the Hamiltonian are 

1, 



E±,„ = (fcn^ + n + -)a;± -\/(A- 2fca;n)2 + Va;2(l + A;n)(l -Fn). (19) 

and the corresponding eigenvectors are 

|-|-, n) = cos6'„|e, n) + sin6'„|g, n -\- 1), (20) 
|-,n) = sm9n\e,n) - cos9n\g,n + 1). (21) 

The expansion coefficients are 



cos^„= , 2^.^(1 + n)(l + fen) 

v/(fin - A„)2 + 492u;2(l + n)(l + kn) 

sin On = ~ ^" = , (23) 

V(f^„ - A„)2 + Va;2(i + + kn) 

in which we have set A„ = A — 2fcnaj and ri„ = A2 + 4g^LU^{l + n){l + kn). 



The energy difference between the levels i?+,n and is ^ A„ + 4:g^oj'^{l + kn){l + n). 
The minimum of the separation occurs when A equals 2knuj and the correspond- 
ing difference is 2g(jj\/ (1 -|- kn){l + n). In Fig. 1 we have plotted the energy 
eigenvalues £'+ and E- as a function of A. The dashed lines represent the 
eigenvalues when g = 0, i.e., £e{g),n- In this case the eigenvalues cross each 
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other as A increases from negative to positive values. The continuous lines rep- 
resent the energy eigenvalues for g = 10~^. The diverging eigenvalue separation 
beyond the minimum separation indicates "level repulsion" in the eigenvalues 
of the dressed atom. The effect of nonzero k is to shift the value of A at which 
the minimum separation or the crossing occurs. If = 0, the minimum as well 
as the crossing occur at A = 0. 



2.2 Evolution of atom-field state 



To understand the dynamics of the atom-field system, we solve for the state of 
the system in interaction picture, where the evolution equation is 

^« = m. (24) 

Here V is the transformed interaction given by 

V = gexp[it{uK+K_ + ^a;,)]{a-K++a+K_)exp[-it{LoK+K_ + ^ct^)]. (25) 

The effect of the transformation on the interaction term is obtained from the 
following results: 

exp{itujK+K)K+ exp{~itLL!K+K^) = K+ eyip{itujKo), (26) 

cxp(-jtj/)cr_, (27) 



itiy itu 
exp( — cr^)CT_ exp( — —a^) 



and their adjoints. Note that the ordering of operators should be maintained in 
the rhs of the first of the results. Using these relations, the interaction picture 
Hamiltonian is written as 

V = g{K+a- exp[it{2(jKo - v)] + adjoint. (28) 

At any time t, let the state of the atom-field system be represented as 



mt)) = E ^e,n{t)\e, n) + Cg,n{t)\g, u) . 



(29) 



n=0 



The coefficients Cg „(t) and Cg.„(t), determined in terms of their initial values 
by the evolution equation Eq.24, are 



exp 



Ce,n{t) — 



. ^flt . ^ Ay^ . 



~2~ 



exp 



iAnt 



Cg,n+l{t) — 



2ig^(n + l)(l + kn) . f il^t 
sm 



2igy^{n+ l)(l + fcn) . 



Ce,„(0) - 

Cg,„+i(0),(30) 

""*^Ce,„(0).(31) 
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The Rabi frequency f2„ is + 45^0;^ (1 + kn){l + n). The dependence of 

fin on n such is that there is a minimum value for the Rabi frequency when n 
satisfies A = 2knuj + g'^iJ^{\ + fc + 2kn){kuj)~^ , provided fc 7^ 0. The variation 
of n„ with respect to n is shown in Fig. 2. In the case of Hjc, the Rabi 
frequency varies hnearly with n and hence there is no minimum. The existence 
of a minimum Rabi frequency has important consequences for the dynamics of 
atomic-inversion, squeezing, photon-statistics, etc. and they are discussed in 
the following sections. 

Let the Rabi frequency attain its minimum for some specific value of n, 
denoted by n. Therefore, we have 

A = A, = 2fca;n+^^^i±^. (32) 

The n dependence of Rabi frequency, for values of n close to n, is obtained by 
Taylor expanding 0„ around n, upto second order. The resultant expression is 

o,..ft, + Srizia!fiv±« (33) 

In Fig. 2 the vahics predicted by the approximate expression for Qn are com- 
pared with those of the exact expression. It is clear that for the chosen values 
of g, k and n, the values match very well for those values of n under the peak 
of the photon number distribution. Quantitatively, the fractional difference is 
less than 3%. 



3 Evolution of atomic inversion 

In the previous section we constructed the complete state of the atom-field sys- 
tem in the dressed atom basis in interaction picture. The state \(t>{t)) in the 
Schrodinger picture is easily obtained by premultiplying the interaction picture 
state function \tp{t)) by exp{i{K-K+ + i'/2)t). In this section, we study the 
temporal evolution of atomic inversion. 

The time-dependent state vector \tp{t)) of the system is determined com- 
pletely once the coefficients Ce,n{t) and Cg^n{t) are known, which, in turn, are 
specified by their initial values. For instance, if the atom is initially in the 
excited state, we have Cg,„(0) = and Ce,ra(0) = {n\il;{0)). The probability 
that the atom is in the excited state, irrespective of the state of the field, is 
I]r=o -^"l^e.nCi)!^ and to be in the ground state is J2^=o Pri\Cg^n{t)\'^ . The 
photon number distribution of the field is P„ . The difference of these two prob- 
abilities is atomic inversion. For the specified initial condition, namely, the atom 
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is initially in the excited state, the atomic inversion W{t) is 



W{t) = l + J2Pn 



Va;2(l+n)(l + fcn) . 

(cos(s2„ij — Ij 



n=0 

and time-dependent part of W{t) is 



fi2 



Writ) = Va;2 ^ P„ii±i!lli±^ cos(f^„f). (35) 

The quantity Wt exhibits rich structure in its evolution. It exhibits collapse 

and revivals when the initial photon distribution is taken to be a Poissonian 
distribution of mean photon number n, which corresponds to the field being in 
a coherent state \\/n). The photon number distribution for the state is P„ = 
exp(-— n)^. This distribution has a single peak and the standard deviation is 
n. Hence, the major contribution to the sum in Eq.35 comes from a few terms 
with n around the peak. With this choice of P„, we have plotted in Fig. 3a the 
evolution of Wt as a function of time. The time required for the first collapse 
and the following revival, denoted by Tc and Tr respectively, can be estimated 
approximately. For the revival to occur, the terms corresponding to those n 
around the peak, should be in phase. Thus, we require Tfi{Cln^i — fin) be equal 
to 27r. The difference of the nearest-neighbour Rabi frequencies is 

2An + A + B 

"n+l - ^'s = ^ • (36) 

The constants A and B are A{k^uj'^ + kg^) and 4(5^ + kg^ — Aktu) respectively. 
The revival time Tr is ^At+A+B ■ 

For the inversion to collapse, the terms in the sum on the rhs of Eq.35 should 
be uncorrelated. Since the width of a Possionian distribution is where the prob- 
ability Pn is appreciable, the condition for collapse is written as Tc{fln_^_^ — 
= 1. If n is large, the expression for Tc is ^J^ . 

The function Wt exhibits rich features when detuning is close to Ac. In Fig. 
3. we have shown the behaviour of Wt for three different values of A. The values 
for the parameters are g = 10~^, k = 10^^, w = 1, n = 30 and the corresponding 
Ac is 0.01606. The values are so chosen that the Rabi frequency attains its 
minimum when n is near n, the average photon number. The evolution of Wt 
with A = .01 < Ac is shown as the top most figure, marked (a) in Fig. 3. The 
figure (6) corresponds to A = Ac and the one marked (c) is for A = .22 > 
Ac. The envelope of Wt when detuning equals Ac, is distinct with structures 
repeating without much distortion. This should be compared with the top and 
bottom figures, which correspond to A 7^ Ac, which exhibit random oscillations 
and do not have neat envelope. 



(34) 
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As noted in Section II, SU(1,1) algebra is realized in terms of K± and Kq 
when k = 1. If we consider resonant interaction (r = 1), then Wt can be 
estimated approximately for the SU(1,1) case. The field is taken to be in a 
coherent state \a) such that |q;| = y/E ^ 1. We use (1 + n)(l + kn) « and 
5^ <C 1 to arrive at f2„ = 2nu}. With these approximations, we arrive at 

Wt = exp(n(cos(cjt) - 1)) cos(n sin(wt)), (37) 

in which we have used (l + n)(l + A;n) ~ f2„^. The magnitude of Wt is negligible 
in this case and so there is no perceptible collapse or revival. This is due to the 

presence of Kerr term and the fact that it dominates over a^a in the Hamiltonian. 
However, we point out that collapses and revivals are present in Wt ii x^^- 



4 Dynamics of field properties 

In the previous section, we studied the dynamics of the two-level atom, in par- 
ticular, the atomic inversion. In the present section, we explore the temporal 
behaviour of field statistics and field amplitudes. As in the previous section, 
the field is initially in a coherent state of complex amplitude a and the atom is 
taken to be in its excited state. With these initial conditions, the probability 
distribution of photons at time t is 



P{n,t) = \Ce,n{t)\^ + \Cg,n{t)\ 
P, 



2 



1 + + ^ cos(n„t) 



(38) 



Coherent states arc the wavepackets whose behaviour is closest to that of a 
classical particle and hence are called classical states. Nevertheless, during 
their evolution in time the states may not be classical. The photon statistics of 
coherent states is Poissonian. Any deviation from this behaviour is characterized 
by Mandel's Q parameter defined as (° ""^"ta/" °^ ■ Using the time-dependent 
probability distribution P(n,t), the expectation values in the expression for Q 
parameter are 

For coherent states, Q is unity. Any value of Q less than unity is nonclassical. 
In Fig. 4 the time evolution of Q parameter for an initially coherent state and 
A = 0.01 < Ac is given. The emergence of nonclassical behaviour {Q < 1) is 
seen. Though not shown in figure, we point out that for A = Ac = 0.016061, 
the statistics does not become sub-Poissonian. When fc = 1 and A = 0, the 
time-dependent part of P{n,t) is of negligible magnitude and so Q does not 
evolve in time. 
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4.1 Squeezing 

We define the field ampfitudes to be 

^W = ^M), (40) 

y^^^^m^^ (41) 

(42) 

These ampfitudes satisfy the commutation relation [X, Y] = i and hence they 
satisfy {6X){6Y) = 1/2. The symbol {6X) stands for the expression ^(^2) - (X)2, 
the variance in X for a given field state. For coherent states of any amplitude a, 
variances in X and Y are the same and equal to 1/2. A state is nonclassical if 
{5X) is less than 1/2, the coherent state value. Using the time-dependent state 
function given in Eq. (30), the variances in the field amplitudes are given by 

{5X f = ^ [1 + (2a^a + a^' + - {a? - {a^? - 2(a^)(a)] , (43) 

and 

{5Yf = \[l + {ra)a - a^^ - a") + {af + {a^f - 2{a^){a)] . (44) 
The expectation values of various operators in these expressions are 

CO 

(a) = Vn+1 [C;„Ce,„+i + C;„C3,„+i] , (45) 

oo 

E 

n=0 



(a') = Y. \/(" +!)(! + kn) [ClnCe,n+2 + C*,„C,,„+2] , (46) 
and 

■oo 

(at a) = ^n[C;„Ce,„ + C;„C,,„]. (47) 

n=0 

The expectation values {a)) and {a)"^) arc the complex conjugates of (a) and 
(a^) respectively. 

The evolution of 5X is shown in Fig. 5. As the field evolves in time, the 
variance in X falls below 0.5 indicating that the quadrature exhibits squeez- 
ing. As a consequence of uncertainty relation, the Y quadrature does not show 
squeezing. However, when we set a = i\/30, the situation is reversed. In this 
case, squeezing is possible in Y and not in X. An interesting feature is that 
the uncertainty in Y with a = i-\/30 is same as in X when a = -s/SO- The 
reason for this is that the field states are of the form X^^g .B„a"|n) wherein 
the coefficients B„ are real for all n. For such superpositions, the uncertainty 
in X for a given a is same as that in Y with a replaced by iQ![25]. 
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4.2 Overlap of initial and time-evolved states 



A quantity of interest is the overlap of the state of the atom-field at time t and 
that at t = 0. With the same initial conditions for the atom- field state as in the 
previous section, the overlap is 

\{mm)}\'=eM-\c^n 

The numerical value of the overlap lies between zero and unity. It is seen from 
Fig. 6 that the overlap becomes zero at longer times. In other words, the 
time-evolved state is almost orthogonal to initial state. For short durations, an 
approximate expression for the decay of overlap is derived by replacing expres- 
sion exp(— n)^ with the Gaussian distribution exp(— '•"gT'' ) '-'^^^ 
n by integration. Further, we set ri„ = fijv + {N — n)il'^ + il'^{N — n)^, with 
the assumption that is close to but not less than n. Here, prime denotes tak- 
ing derivative with respect to n and the suffix represents the value of n where 
the derivative is evaluated. With these approximations we get, after neglecting 
oscillatory terms, 

E°° (n-n)2 Mnts 

n=0 

OC 



Similar expression can be derived for summation with sin(^^). The above 

expression indeed predicts that the overlap function decays with time. When 
N = n, the first derivative of f2„ vanishes and the exponential term in the 
envelope is absent. Consequently the decay is slower. However, if the photon 
number distribution of the field is very broad, it is incorrect to truncate the 
Taylor series and the expression in Eq. (49) is invalid. 

5 Summary 

The Hamiltonian that contains the usual and intensity-dependent (including 
Kerr term) JC models as limiting cases has been constructed. The nondissipa- 
tive dynamics dictated by the generalized Hamiltonian H is completely solvable. 
The eigenvalues of the Hamiltonian, with Kerr term and intensity-dependent 
coupling, exhibits level-repulsion. In the case of nonvanishing Kerr term, the 
Rabi frequency, considered as a function of the photon number n, attains a 



E 

n=0 



,|2n 



n! 



cos(— ) + 



(48) 



Re 



1 + N' 
X exp 



o^n 4 
dn 4 



(49) 
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minimum. The dynamical behaviour of atomic-inversion, when detuning is so 
chosen to have the minimum Rabi frequency, exhibits superstructures, which is 
absent in the usual JC model. The overlap of the initial coherent state and the 
time-evolved state decays with time. In the language of inner product, the ini- 
tial coherent state is almost orthogonal to the time-evolved state. The results, 
except the occurrence of superstructures in atomic inversion, a consequence of 
nonvanishing Kerr term, go over to those of the JC model, in the limit of A; — » 0. 

We note that the formal equivalence between Hjc and H is obtained by 
identifying with and K- with d. This, in conjunction with the fact both 
the sets of operators {d^,d, 7} and {K±,Ko} are closed under commutation, 
implies that the expansion coefficients Ce/g,,i(t) for the evolution governed by 
H are obtained from the corresponding expressions for the usual JC model by 
replacing A by A„ and g by gV^ + kn. Hence, all those physical quantities, 
like the atomic inversion, quadrature fluctuations, etc., computed in terms of 
the expansion coefficients are derivable from those of the JC model. 
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Fig. 1 Dependence of eigenvalues and E- on detuning A. The continuous 
curve corresponds to 5 = 0.1 and the dashed curve corresponds to g = 0. Here 
k = .1 and n = 30. The dashed curves with positive and negative slopes corre- 
spond respectively to fe,n and £g,n- Lower part of the figure is for n = 1 and 
upper part for n = 2. 

Fig. 2 Variation of ri„ with n. We have set g = 10~^, k = 10"'', n = 30 
and A = 0.016061. The approximate and actual Rabi frequencies are compared 
in the upper figure. Dotted line corresponds to the approximate expression in 
Eq.33 and continuous curve corresponds to exact expression. Values of k, g, n 
and A are 10~^, 10~^, 30 and 0.016061 respectively. The bottom curve shows 
the photon number distribution for the coherent state \a = \/30). 

Fig. 3 Time-dependent part of atomic inversion. Case (a) corresponds to 
A = 0.01 < Ac. Case (b) is for A = A^ = 0.016061 and Case (c) refers to 
A = 0.022 > Ac. Here Case (b) shows 1 + Wt and Case (c) shows 2 + Wt- 

Fig. 4 Time variation of {SX)"^, as the atom-field system evolves. The evolu- 
tion is shown for three values of detuning, 0.01 (continuous), 0.016061 (dotted) 
and 0.02 (dashed). Instants of {SX)"^ less than 0.5 correspond to squeezing in 
X quadrature. 

Fig. 5 Mandel's Q parameter as a function of time. Instants of Q < 1 cor- 
respond to sub-Poissonian statistics. Values of g, k and n are same as in Fig. 2 
and A = 0.01. 



Fig. 6 Overlap of initial and time-evolved field states. Y-axis corresponds to 

|(^(0)|^(t))p. The envelope of the overlap function decays with time implying 
that the initial state is almost orthogonal to the evolved state. The values of 
the parameters are same as in Fig. 2 and the detuning A is equal to 0.016061. 
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